5 Forced oscillators and limit cycles 



5.1 General remarks 



How may we describe a forced oscillator? 
The linear equation 



e + ^e + u'^o = 



(3) 



is in general inadequate. Why? 

Linearity =^ if 6{t) is a solution, then so is a6{t), a real. This is incompatible 
with bounded oscillations (i.e., ^max < tt). 

We therefore introduce an equation with 

• a nonlinear ity; and 

• an energy source that compensates viscous damping. 

5.2 Van der Pol equation 

Consider a damping coefficient 7(^) such that 



7(^) > 



for 16*1 large 



7(^) < for 1^1 smah 



Express this in terms of 9^: 




where 70 > and ^0 is some reference amplitude. 



Now, obviously. 



7 > 
7 < 



for > el 
for 0^ < el 
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Substituting 7 into (3) , we get 




This equation is known as the van der Pol equation. It was introduced in the 
1920's as a model of nonUnear electric circuits used in the first radios. 

In van der Pol's (vaccum tube) circuits, 

• high current =^ positive (ordinary) resistance; and 

• low current =^ negative resistance. 

The basic behavior: large oscillations decay and small oscillations grow. 

We shall examine this system in some detail. First, we write it in non- 
dimensional form. 

We define new units of time and amplitude: 

• unit of time = l/u 

• unit of amplitude = ^o- 

We transform 



e 



t 



t'/u 



where 6' and t' are non-dimensional. 



Substituting above, we obtain 




Divide by u'^do: 




dt'2 
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Now define the dimensionless control parameter 



Finally, drop primes to obtain 



What can we say about the phase portraits? 

• When the amplitude of oscillations is small (6'niax < 1), we have 

^(^max ~ 1) < =^ negative damping 
Thus trajectories spiral outward: 



(4) 




• But when the amplitude of oscillations is large (^max > 1); 

■^(^max ~ 1) > ^ positive damping 
The trajectories spiral inward: 

e 



e 
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Intuitively, we expect a closed trajectory between these two extreme cases: 



This closed trajectory is called a limit cycle. 

For e > 0, the limit cycle is an attractor (and is stable). 

This is a new kind of attractor. Instead of representing a single fixed point, 
it represents stable oscillations. 

Examples of such stable oscilations abound in nature: heartbeats (see Fig- 
ure from Glass); circadian (daily) cycles in body temperature, etc. Small 
perturbations always return to the standard cycle. 

What can we say about the limit cycle of the van der Pol equation? 

With the help of various theorems (see Strogatz, Ch. 7) one can prove the 
existence and stability of the limit cycle. 

We may, however, make substantial progress with a simple energy balance 
argument. 

5.3 Energy balance for small e 

Let £ — ^ 0, and take 9 small. Using our previous expression for energy in the 
pendulum, the non-dimensional energy is 
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The time variation of energy is 

From the van der Pol equation (4), we have 

e = -£{e^ - 1)9 - 9. 
Substituting this into the expression for dE/dt, we obtain 



= e9\l - 9') - 99 + 99 (5) 

= e9\l-e^) (6) 



Now define the average of a function f[t) over one period of the oscillation: 

Then the average energy variation over one period is 

dE 1 d^ 



dt 27r Jt^ dt 



■dt. 



Substituting equation (6) for dE/dt, we obtain 



dE 



= £9'^ - e9'^9'^. 



dt 

In steady state, the production of energy, 69"^, is exactly compensated by the 
dissipation of energy, s9'^9'^. Thus 



or 



^2 = ^2^2 _ 



Now consider the limit e ^ (from above). 
We know the approximate solution: 

9{t) = /osint. 
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i.e., simple sinusoidal oscillation of unknown amplitude p 
We proceed to calculate p from the energy balance. 
The average rate of energy production is 

p^ cos^tdt = -p 

The average rate of energy dissipation is 

p^ sin^ t cos^ tdt = -p 



1 /'to+27r 1 



2ti 



'to 



1 pto + 2TT I 

0202 ~ _ / p4 ^ ^Qg2 ^ _^„4 



27r 



'to 



The energy balance argument gives 



Therefore 



1 2 1 4 

2^ = 8^ • 



p = 2. 



We thus find that, independent of e = 70/cj, we have the following approxi- 
mate solution for £ <C 1: 

e{t) ~ 2sint. 

That is, we have a limit cycle with an amplitude of 2 dimensionless units. 
Graphically, 



e 



e 



CD 




Further work (see, e.g., Strogatz) shows that this limit cycle is stable. 
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5.4 Limit cycle for e large 



The case of e large requires a different analysis. We follow the argument 
given in Strogatz (p. 212). 

First, we introduce an unconventional set of phase plane variables (not x = 
y,y = . . .). That is, the phase plane coordinates will not be 6 and 6. 



Recall the van der Pol equation (4), but write in terms oi x = 6: 

X + e{x — l)x + X = 0. 

Notice that 



(7) 



X + ex(x'^ — 1) = — 
^ ^ dt 



^ 3 

X + £ -X — X 



Let 



and 



F{x) = ^x^ — X 



w = X -\- eF{x). 
Then, using (8) and (9), we have 

w = X -\- ex{x'^ — 1). 

Substituting the van der Pol equation (7), this gives 

w = —X 

Now rearrange equation (9) to obtain 

X = w — eF[x) 



(8) 
(9) 



(10) 



(11) 



We have thus parameterized the system by x and w. However we make one 
more change of variable. Write 



Then (10) and (11) become 



y = w/e. 



X = s[y — F{x)] 
1 

y = — X 

e 



(12) 
(13) 
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Now consider a trajectory in the x-y plane. 

First, draw the nullcline for x, that is, the curve showing where x = 0. This 
is the cubic curve y = F{x). 




Now imagine a trajectory starting not too close to y = F{x), i.e.. suppose 

y-F{x)^l. 
Then from the equations of motion (12) and (13), 

X ^ e ':$> 1 

y ~ 1/s^l assuming x ~ 1. 

Thus the horizontal velocity is large and the vertical velocity is small. 
^ trajectories move horizontally. 

The y-nuUcline shows that the vertical velocity vanishes for x = 0.) 

Eventually the trajectory is so close to y = F[x) such that 

y - F{x) ~ \ 

implying that 

X '-^ y —. 

£ 

Thus the trajectory crosses the nullcline (vertically, since i = on the null- 
cline) . 

Then x changes sign, we still have ± ~ ^ ~ 1/e, and the trajectories crawl 
slowly along the nullcline. 
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What happens at the knee (the minimum of F{x))7 

The trajectories jump sideways again, as may be inferred from the symmetry 
-X, -y. 

The trajectory closes to form the Umit cycle. 

Summary: The dynamics has two widely separated time scales: 

• The crawls: At ~ e: {x ^ 1/e) 

• The jumps: At ^ 1/e {x ^ e) 

A time series of x(t) = 0(t) shows a classic relaxation oscillation: 



Relaxation oscillations are periodic processes with two time scales: a slow 
buildup is followed by a fast discharge. 

Examples include 

• stick-slip friction (earthquakes, avalanches, bowed violin strings, etc.) 

• nerve cells, heart beats ( large literature in mathematical biology...) 

5.5 A final note 

Limit cycles exist only in nonlinear systems. Why? 
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A linear system x = Ax can have closed periodic orbits, but not an isolated 
orbit. 

That is, linearity requires that if x{t) is a solution, so is ax{t), a 0. 

Thus the amplitude of a periodic cycle in a linear system depends on the 
initial conditions. 

The amplitude of a limit cycle, however, is independent of the initial condi- 
tions. 
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